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ABSTRACT 

The explosion of a carbon-oxygen white dwarf as a Type la supernova is known to be sensitive to the manner 
in which the burning is ignited. Studies of the pre-supernova evolution suggest asymmetric, off-center ignition, 
and here we explore its consequences in two- and three-dimensional simulations. Compared with centrally 
ignited models, one-sided ignitions initially burn less and release less energy. For the distributions of ignition 
points studied, ignition within two hemispheres typically leads to the unbinding of the white dwarf, while 
ignition within a small fraction of one hemisphere does not. We also examine the spreading of the blast over the 
surface of the white dwarf that occurs as the first plumes of burning erupt from the star. In particular, our studies 
test whether the collision of strong compressional waves can trigger a detonation on the far side of the star as 
has been suggested by Plewa et al. (2004). The maximum temperature reached in these collisions is sensitive 
to how much burning and expansion has already gone on, and to the dimensionality of the calculation. Though 
detonations are sometimes observed in 2D models, none ever happens in the corresponding 3D calculations. 
Collisions between the expansion fronts of multiple bubbles also seem, in the usual case, unable to ignite 
a detonation. "Gravitationally confined detonation" is therefore not a robust mechanism for the explosion. 
Detonation may still be possible in these models however, either following a pulsation or by spontaneous 
detonation if the turbulent energy is high enough. 

Subject headings: Stars: supernovae: general — hydrodynamics - instabilities — turbulence — methods: 
numerical 



1. INTRODUCTION 

In the currently favored model for Type la supernovae 
(SNe la), a carbon-oxygen white dwarf (WD) grows to al- 
most the Chandrasekhar mass, then explodes due to a ther- 
monuclear instability (e.g., [Hillebrandt & Niemever 2000). 
While the modeling of the explosion itself has reached a 
high level of sophistication, with multi-dimensional stud- 
ies being routinely carried out by several groups (e.^ 



i lGamezo et al. |200l 'Plewa et a lJ2004l 
I2005at iGarci'a-Senz & Brave 



Reinecke et al .'20021: 
^ooke & Hillebrandi 12004 
2005t iRopke et a l. 2006a'), the initial conditions of this pro- 
cess remain largely unknown. This is unfortunate since the 
geometry of the fla me ignition has a large effect on the 
explosion s trength ( Niemeyer et al. 'H^. Ch^ne etai]l2(M 
lR5pke et al. 2006b; Schmidt & Niemeyer 2006). 

Runaway commences once the WD has accreted sufficient 
matter from a binary companion to approach a central density 
~ 3 X lO' g cm""* where plasma neutrino losses are exceeded 
by energy generation from a highly screened carbon fusion 
reaction. The stage for flame ignition is set by a century of 
convective carbon burning in the progenitor WD. It remains 
unclear, however, if the first sparks to develop a nearly discon- 
tinuous temperature gradient on their perimeters ( the "flame") 
are c oncentrated in the center of the star ( Hoflich & SteinI 
l2002h or spread around by the convective flow in which 
they are embedded. Since this convective flow may have a 
dipole character, one natural possibility is lopsided ignition 
displaced somewhat off-center (e.g., W oosley et a l. 2004). 

It is known that, once born, the flame experiences an 
extended period of subsonic propagation — a "deflagration" 
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JNomoto et alJIlQT^ . Prompt detonation is excluded on the 
grounds that it would produce spectroscopy, nucleosynthe- 
sis, and a light curve very different from observations (see 
iFilippenko 1997, for a review). It would also require a de- 
gree of isot hermal! ty in the core that would be very difficult 
to achieve llWooslevlll990l) . The deflagration poses a com- 
putational challenge since the ashes of the burning are buoy- 
ant, and that leads to instabilities and turbulence that can only 
be followed with any accuracy in a multi-dimensional cal- 
culation. The difficulty is compounded by the large range 
of spatial scales — sub-millimeter for the flame width and 
Kolmogorov scale to 2,000 km for the WD — and the high 
Reynolds number. Re ~ 10'"*. 

Calculations of the same WD differ in outcome because 
of the assumptions about ignition, various techniques used 
to treat flame instabilities, and turbulence. It is not practi- 
cal to resolve both the flame and the star, so full-star models, 
such as the ones presented here, rely upon an effective flame 
model and a subgrid scale model for the turbulence. Quali- 
tatively, the flame advances radially at a speed given by the 
flotation of the largest plumes, but the lateral spreading of 
each plume and the overall efficiency of the explosion can 
vary, depending upon the way turbulence is handled, and on 
the dimensionality and resolution of the model. Results of dif- 
ferent approaches have been published in a variety of studies 
with the general conclusio n that a pu re deflagration can give a 
viable explosion (Reineck e et alJl20 02b: Gamezo et al. 20021 
iRopke & Hillebrandt 2005a'), not too different from what is 
observed (Travaglio et al. 2004; Blinnikov et al. 2006). It re- 
mains controversial, however, if these models can give light 
curves as bright as some observations indicate, or can ex- 
plain all of the s pectroscopic features dGamezo et alJl2003l 
lKozmaetal.l2(i05). 

Moreover, these successful models all have in common the 
assumption of nearly isotropic central ignition. It may be 
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that nature provides only anisotropic ignition conditions, or 
it may be that the the observational constraints on a pure de- 
flagration will ultimately prove too stringent. In these cases, a 
transition to det onation may need to occur (Plewa et al. 2004; 
iLivne et alf 2005>. The idea of a delayed detonation has been 
aroun d for some time (Khokhlov 1991; Wooslev & Weaver 
[1994), but the physics of that transition, if it happens, is still 
uncertain ( Niemever & WooslevMl997;iNiemever.il99^ . Re- 
cent two-dimensional calculations have suggested that burned 
material may quickly ascend to the surface of the still gravi- 
tationally bound star, sweep around it and, by collision and 
compression on the opposite side, tr igger a detonation in the 
unburned material iPlewa et alJl2004i) . This is called by its 
proponents "gravitationally confined detonation", or GCD. 

Here, we follow the evolution of one-sided ignitions for a 
range of assumptions regarding the initial conditions in two- 
dimensional (2D) and three-dimensional (3D) simulations. 
The results are sensitive to the fuel consumption and energy 
release during the rising stage of the plume. Therefore a cor- 
rect description of the turbulent deflagration flame is as crucial 
here as in other models. We describe our approach in Sect.|3l 
Before that, we give a brief motivation o f the ignition scenar- 
ios explored here. Although lRopke et al.l(.2006b.) showed that 
ignition conditions cannot be explored reliably in 2D simu- 
lations, we can and do use surveys in cylindrical symmetry 
to get a feeling for the parameter range to be explored, as 
well as the dependence of the results on numerical resolution 
(Sect. 13. Full-star 3D simulations are presented in Sect. |5j 
and the consequences regarding the possibility of triggering a 
detonation are discussed in Sect. |S] Conclusions are drawn in 
Sect.lll 

2. OFF-CENTER IGNITION 

Several studies now suggest ignition with an offset 
from the center of the WD of order 100... 200 km 
jOarcia-Senz & Wooslev |1995|; jpooslev et al. 20041 
IWunsch & Wooslev 2004; Kuhlen et al. 2006; lapichino et al. 
12006). It should be acknowledged that none of these studies 
has yet followed the actual transition from a high temperature 
fluctuation to a flame in a self-consistent way, including the 
possibility the perturbation is disrupted by turbulence, and all 
fall far short of the actual Reynolds number in the star Also, 
while arguments based upon a probability density function 
can offer some guidance as to whether a particular tempera- 
ture is likely to be realized, they cannot, by themselves, say 
whether the high temperature happens in a contiguous region, 
or in many disparate points, or even over some interval of 
time. 

Still, the calculations of lKuhlenet al. (2006) do suggest that 
ignition is unlikely to occur as a single spherical bubble either 
at the center or off-center. Rather the distribution of high tem- 
perature may look more like a "teardrop", spreading as it goes 
out to a large opening angle. We thus explore here a variety 
of initial conditions ranging from nearly spherical bubbles far 
off center, to ignition in multiple points forming complex con- 
figurations. The models are named according to their dimen- 
sionality and the ignition characteristics where "B" stands for 
a single bubble (spherical as much as Cartesian coordinates al- 
low) with an attached number indicating its radius in kilome- 
ters, "P" is a highly perturbed bubble, "T" stands for a teardrop 
distribution of multiple bubbles followed by a "1" to indicate 
a strictly one-sided ignition and by a "2" when the ignition re- 
gion overshoots to the opposite side. Where appropriate, "d" 
followed by a number gives the distance of the ignition center 



from the middle of the WD in the bubble case (and the maxi- 
mum extent of ignition measured from the WD's center in the 
teardrop case). For example. Model 3B50d200 is a 3D simu- 
lation with an ignition in form of a spherical bubble of radius 
50 km centered 200 km from the middle of the star There 
may also be variations on these names based upon resolution 
(a,b,c, etc). 

3. THE ASTROPHYSICAL AND NUMERICAL MODEL 

The implementation of the deflagration SN la model fol- 
lows the detailed descriptions given by Reinecke et al. ( 199^ 
2002a); Ropke (2005) and Schmidt etal. (2006b). The hy- 
drodynamics is modeled via the piecewise parabolic method 
(Colella & Woodward 1984) in the Prometheus implementa- 
tion ( Frvxell et a l. 1989), in combination with a WD matter 
equation of state incorporating an electron gas relativistic and 
degenerate to variable degrees, a Boltzmann gas of nuclei, 
photons and electron-positron pairs. 

Since the results are sensitive to small perturbations of the 
initial flame configuration, the directional splitting in the hy- 
drodynamics solver, as well as discretization errors on the 
Cartesian grid may provide seeds for developing instabilities. 
This is not necessarily unphysical. In a realistic SN la explo- 
sion, the background is not expected to be smooth, nor will 
the initial flame shape be perfectly regular. But, as discussed 
below, one sometimes has little control over these effects. 

Another crucial aspect of the modeling is the prescription 
for flame propagation. We strive to implement here a consis- 
tent model for burning in the flamelet regime following stan- 
dard theories of turbulent combustion — such as they are. The 
fundamental assumption is the establishment of a turbulent 
cascade from large-scale eddies produced by shear instabil- 
ities at the interfaces of burning bubbles down to the Kol- 
mogorov scale, where viscous effects dissipate the turbulent 
kinetic energy. The flame interacts with turbulence down to 
the Gibson scale, where the laminar burning speed s\ becomes 
comparable to the turbulent velocity fluctuations u'. By defi- 
nition in the flamelet regime, this Gibson scale is large com- 
pared to the width of the flame — a condition that holds for 
most parts of the supernova explosion until quite low densi- 
ties are reached ( for an approach to modelin g stages beyond 
this regime see Ropke & Hillebrandt 2005b). Thus the inter- 
nal flame structure remains unaffected by turbulent eddies. 
Damkohler ( 194^ first pointed out that the turbulent flame 
front in these circumstances should propagate with an effec- 
tive velocity that is proportional to the turbulent velocity fluc- 
tuations and independent of the burning microphysics. 

Following this concept, the flame is modeled as a sharp 
interface separating the fuel from the ash es. Its propaga- 
tion is followed in a level set approach C Osher & Sethianl 
1988), where the flame front is associated with the zero level 
set of a scalar, G, defined to be a signed distance function 
away from the flame. This front is advanced as described by 
Reinecke et al. (1999). 

The speed of flame propagation is a function of the turbu- 
lent velocity fluctuations on the scale of the computational 
grid cells. Its value is determined from a subgrid-scale tur- 
bulence model. In the 2D simulations, this model is imple- 
mented according to Niemever & Hillebrandt ( 1995), and the 
turbulent flame speed, s^, is set equal to the turbulent velocity 
fluctuations, u' = sjlk^^, where k^gs denotes the subgrid-scale 
turbulent specific kinetic energy. The 3D simulations employ 
an improved subgrid-scale model ( Schmidt et al. .2006a h) . 
That model, based upon localized closures for the terms of 
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the balance equation of turbulent subgrid-scale energy, does 
not need to assume a specific scaling behavior of the turbulent 
cascade, nor isotropy of the turbulence. The flame propaga- 
tion speed is implemented as 



with Ct = 4/3 and (7si>s denoting the subgrid-s cale turbulence 
velocity ('Pocheaul994'; 'Schmidt et al.l2006b »). 

Nuclear reactions are implemented in a simple way de- 
scribed by Rei necke et alJ (l20()2all) . Only four species, '^C, 
'^O, ^'^Mg, ^^Ni, and a-particles, are taken into account. Ma- 
terial traversed by the flame front is converted to a com- 
position represented by a temperature-dependent mixture of 
^^Ni and a-particles in nuclear statistical equilibrium, or, 
at fuel densities below 5.25 x 10^ gem"-', to intermediate 
mass elements represented by ^''Mg. At fuel densities below 
10^ g cm"^ nuclear reactions are assumed to cease. 

In all simulations the WD was set up cold and isothermal 
with a temperature of T = 5 x 10^' K and a central density of 
Pc = 2.9 X lO' g cm"^ composed of equal parts of carbon and 
oxygen. 

The discretization on the computational grid follows the 
strategy of two nested moving grids suggested bvl Ropke et al.l 
(|2006b), where a fine-resolved uniform inner grid contains 
the flame while an outer grid with exponentially growing grid 
cells accommodates the WD star and follows its expansion. 
Due to flame propagation inside the WD, it is possible to sub- 
sequently collect adjacent grid cells of the outer grid into the 
uniform part, as soon as the cell sizes match thereby optimiz- 
ing the resolution for the given number of computational grid 
cells. 

4. EXPLORING THE POSSIBILITIES: 2D SIMULATIONS 

The parameter space and the dependence of the results on 
numerical resolution were first explored assuming cylindrical 
symmetry. The results obtained in 2D should not be used to 
draw quantitative conclusions, but are a numerically inexpen- 
sive way to explore a broad range of possibilities. 

In the cylindrical (r, z)-setups, the entire WD was accom- 
modated on the grid, from its center to its radius assuming 
rotational symmetry about the z-axis (see top left panel of 
Fig.Q. 

4.1. Single-bubble ignition 

As a first numerical experiment, the flame was initiated 
in the simplest conceivable configuration: a single spherical 
bubble ignited somewhere on the z-axis. Even this simple 
configuration has three parameters that potentially impact the 
evolution of the explosion. One is the resolution of the flame 
and the WD star. Two additional parameters — the displace- 
ment of the igniting bubble from the center of the WD and its 
radius — are more physical in nature. As an illustrative model, 
consider one where ignition took place 200 km off'-center in a 
spherical bubble of radius 50 km (cf. Fig.Q. The resolution 
for this example study (Model 2B50d200c, Tabled was 256 
grid cells in r-direction and 512 grid cells in z-direction, and 
the initial setup is shown in the top left panel of Fig. ^ In 
the top middle panel, a close-up on the flame illustrates the 
evolution from the initial spherical shape to a more irregular 
toroidal structure. This is a consequence of the buoyancy- 
induced flotation of the bubble acting in combination with 
self-propagation of the flame due to burning. The flame floats 



towards the surface of the star, thereby being subject to con- 
siderable lateral spread (cf. snapshots at f = 0.5 s and f = 1.0 s 
in Fig.0. At f ~ 2 s, burning has ceased since the fuel density 
ahead of the flame has dropped below 10^ g cm"'. The burned 
material sweeps around the star since it is still gravitationally 
bound. Although it extends to rather large radii, only about 
0.49 Mq is located outside a radius of 2.5 x 10^ cm at f = 2 s. 
Only ~ 0.1 Mq of the WD is burned, inadequate to unbind 
the star In the outer layers of the WD opposite to the flame 
ignition, burned material collides subsonically with a com- 
pressional front moving ahead of the actual ash. In the sim- 
ulation, the evolution continues by an expansion of the outer 
layers while the central parts of the WD contract such that the 
inner parts of the flame eventually reach fuel densities above 
the threshold of 10^ g cm""*. In our flame description burning 
resumes at this point, and the newly processed material is ex- 
pelled into the ash region (cf. Fig. ^ snapshot at f - 6.0 s). 
This, however, is sensitive to the way burning is implemented 
in the code (either completely "on" for densities above 10^ g 
cm"^, or "off" for lower densities) and may not be a realistic 
occurrence. 

As pointed out by Plewa et al. (2004), the collision of the 
gravitationally bound material sweeping around the surface 
of the WD marks an interesting point in the evolution. To 
initiate a spontaneous detonation, high density and tempera- 
ture are both necessary in the unburned material at the focus 
of the collision. The cylindrical-symmetric setup forces the 
collision of burned material to take place at the negative part 
of the z-axis. Figure |2l shows the temperature of the material 
in the model described above. The snapshot was taken at the 
instant of peak temperature in the compressed fuel in the col- 
lision region. The peak temperature there was 2.22 x 10^ K 
at a density of 1.41 x 10'' gem"-'. As will be discussed, the 
conditions reached in the collision depend on the parameters 
of the setup as discussed. 

4.1.1. A note on numerical convergence 

Two sets of simulations with spherical initial flames were 
carried out varying the number of cells in the computational 
grid. One had an initial bubble radius of 50 km, and the other, 
a bubble radius of 25 km. For each simulation, the maximum 
temperature in the unburned material in the collision region 
was determined, along with the nuclear energy release prior 
to reaching this temperature. The results are given in Tabled 
There is significant scatter in the critical collision temperature 
with values deviating by up to 43%. A similar variation is 
seen in the energy of the burning. Although a clear trend is not 
apparent in Tabled it seems that better numerical resolution 
leads to less energy release in the burning, less expansion, and 
hence to a stronger collision in the unburned material. The 
scatter is larger when starting with bubbles of smaller radius. 

Previous studies suggest convergent results should be ob- 
tained with a resoluti on of about 256 cells per dimension in 
one octant jReinecke et al J2002at iRopke 2005). Here the sit- 
uation is different. While previous resolution studies were 
carried out on the basis of an initial flame setup with well 
specified perturbations imposed on it, all seeds for growing 
nonlinear instabilities in the spherical bubble setup applied 
here are introduced by numerical artifacts, such as discretiza- 
tion errors and noise. This setup is not a well-posed numeri- 
cal problem because the evolution of nonlinear features in the 
flame structure is expected to be strongly resolution depen- 
dent and thus numerical convergence is problematic. Given 
that the flame evolution is dominated by strongly nonlinear 
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Fig. 1 . — Evolution of a 2D explosion simulation (Model 2B50d200c) ignited in a single bubble 200 km oflF-center. The top-middle panel illustrates the flame 
front evolution in the interval ; = [0, 0.3] s. Each contour con'esponds to a time step of 0.05 s. In all other snapshots, the cyan isosurface con'esponds to the zero 
level set of G which is associated to the flame front in early stages of the evolution and indicates the interface between fuel and ashes once burning has ceased. 



effects, the variation in Table 1 is no big surprise and is illus- 
trative of the uncertainty in our results. We emphasize though 
that this is not due to the numerical methods applied here, but 
due to the variable (and artificial) initial setup. 

Viewed this way, the results of the resolution study can be 
understood in a straightforward manner It should be noted 
that the resolution affects the answer in two ways. On the 
one hand, it affects the temperature in the collision because 
higher resolution smears out the hot spot less. But, on the 
other hand, it also affects the propagation of the flame as it 
moves outwards in the star, as illustrated by the variable ex- 
plosion energy. Higher resolution decreases the discretization 
errors and therefore reduces the seeds for the growth of non- 



linear perturbations. Therefore the flame develops less sur- 
face, less material is consumed and the lower energy release 
leads to a weaker expansion of the star The material sweeping 
around at the surface is stronger gravitationally bound and the 
clash is more vigorous. Smaller initial flame bubbles are more 
sensitive to the numerical resolution. This is understandable 
since here bubble flotation is slower and nonlinear features in- 
voked by discretization errors have more time to develop on 
the way to the surface of the WD. 

An imaginary ideal situation with no discretization errors 
would suppress the nonlinear growth due to instabilities. But 
since the WD star is expected to be perturbed by pre-ignition 
convection and the flame is likely to ignite in multiple spots 
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TABLE 1 

2D RESOLUTION STUDY FOR IGNITION 200 km OFF-CENTER IN A BUBBLE OF 50 km RADIUS. DaSHES MARK CASES WHERE THE DENSITY AT WHICH TO MEASURE THE TEMPERATURE 

IS NOT REACHED IN THE COLLISION. 



Model 


bubble 


resolution 


r„,a.K at coll. 


Eauc at coll. 


rmax(p>3xl06gcm-3) 


rmax(p> lxlO^gcm-3) 


surface 




radius [km] 




[10' K] 


[10^" erg] 


at coll. [10' K] 


at coll. [10' K] 


detonation? 


2B50d200a 


50 


128 X 256 


2.61 


1.14 


1.54 




no 


2B50d200b 


50 


192 X 384 


2.92 


0.97 


2.60 




yes 


2B50d200c 


50 


256x512 


2.22 


1.46 


1.28 




no 


2B50d200d 


50 


384 X 768 


2.53 


1.44 


0.959 




no 


2B50d200e 


50 


512 X 1024 


2.29 


1.39 


0.954 




no 


2B25d200a 


25 


128 X 256 


2.40 


1.33 


2.08 




no 


2B25d200a 


25 


192 X 384 


1.97 


1.47 


0.224 




no 


2B25d200a 


25 


256x512 


2.60 


1.09 


2.32 




yes 


2B25d200a 


25 


384 X 768 


3.03 


0.72 


3.03 


2.95 


yes 


2B25d200a 


25 


512 X 1024 


3.83 


0.82 


3.83 


3.80 


yes 




2 0-1 [)°4.0-1CrS.0-l[fB.0-10^. 0-1 0'l. 2-10' 

X [cm] 

Fig. 2. — Temperature distribution Model 2B50d200c. The white contour 
indicates the zero-level set of G indicating the interface between fuel and 
ashes. Apart from the ash regions enclosed by this contour, significantly 
increased temperatures are found in the region compressed by the collision 
(marked by the white arrow). 

or an irregular shape, this seems far from reality. 

Besides the energy release due to burning, the impact of the 
colliding material will also be sensitive to the morphology of 
the ash that is driving the unburned material like a piston. A 
well-defined, large leading edge of the colliding ash regions 
should result in a better focus than a multitude of leading fea- 
tures. Since the flame morphology is determined by nonlinear 
effects and instabilities, discretization errors in diff'erent reso- 
lutions amplify the scatter in the results. 

4.1.2. Bubble displacement 

The first physical parameter — the displacement of the ini- 
tial flame bubble from the center — turned out to have sub- 
stantial influence on the densities and temperatures reached 
in the collision. Values from simulations starting with an ini- 
tial flame bubble of radius 25 km and varying distances from 
the center of the WD are given in Table |2] In order to min- 
imize the scatter due to diff'erences in resolution, the initial 
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Fig. 3. — Maximum fuel temperature reached in the region of collision of 
the material sweeping around the WD as a function of the displacement of 
the initial flame bubble from the center of the WD (according to Tablel2l. 



grid spacing of the 256 x 512 cells setup was held fixed for 
all simulations. Note, however, that since two nested grids 
are used to follow the expansion of the WD and the flame 
propagation, the resolution evolves according to the energy 
released in the burning, which is diff'erent in the various sim- 
ulations. Nevertheless, a trend of increasing collision temper- 
atures with larger initial flame displacements from the center 
of the WD is clearly visible (cf. Fig. |3}- At large displace- 
ments the temperature increases less indicating a saturation 
of the eff'ect. A possible explanation is that although here 
the expansion of the star prior to breakout is decreased, the 
amount of ashes expelled from the surface also decreases due 
to less burning taking place. Therefore the momentum of 
the colliding ash regions is smaller and the reduced impact 
leads to lower compression temperatures. Flames born closer 
to the center burn more material and cause more expansion 
(cf. Tabled}- Displacing the initial flame bubble from 50km 
to 600 km off'-center decreases the nuclear energy release for 
94% and increases the maximum temperature reached in the 
collision by 250%. Such a large displacement as 600 km is 
not realistic, but the consequences may be the same as for a 
bubble ignited closer in, but with a less efficient prescription 
for burning on the way out. 
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TABLE 2 

Collision parameters in dependence of initial bubble displacement (2D simulations; initial bubble radius: 25 km). Dashes mark cases where the density at 

WHICH TO measure THE TEMPERATURE IS NOT REACHED IN THE COLLISION. 



Model 


distance from 
center [km] 


r^ax at coll. 
[10' K] 


£nuc at coll. 
[lO^^erg] 


7'max(p>3xl0«gcm-3) 
at coll. [10' K] 


7'max(p> lxlO'gcm-3) 
at coll. [10' K] 


surface 
detonation? 


2B50d50 


50 


1.46 


2.31 


0.809 




no 


2B50dl00 


100 


2.00 


1.93 






no 


2B25d200 


200 


2.81 


1.07 


2.13 




no 


2B25d300 


300 


3.45 


0.76 


3.45 


0.438 


yes 


2B25d400 


400 


3.91 


0.46 


3.91 


3.91 


yes 


2B25d500 


500 


4.54 


0.28 


4.54 


3.91 


yes 


2B25d600 


600 


4.61 


0.12 


4.61 


4.61 


yes 
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Fig. 4. — Different initial flame bubble morphologies. Left: perfect sphere 
(dashed) vs. irregular shape (solid); right: teardrop- shaped initial flames on 
one side of the WD's center (dashed) and overshooting to the opposite side 
of it (sohd). 

4.1.3. Bubble morphology 

The initial morphology of the flame also affects the strength 
of the collision. Any divergence from a perfectly spherical 
shape has a similar effect to varying the resolution. In both 
cases, seeds for the developing instabilities are imposed — 
with changing resolution due to discretization errors, and for 
more complex initial bubble shapes, explicitly in a controlled 
way. 

To demonstrate this, two simulations were carried out on 
a well resolved (384 x 768 cells) computational grid. One 
model (2B50dl00, see Table |3j was ignited in a — within dis- 
cretization error — ^perfectly spherical bubble of radius 50 km 
at a distance of 100 km from the center. In the second simula- 
tion (2P50dl00), the initial flame was composed of 160 par- 
tially overlapping small bubbles of radius 3 km placed within 
a sphere of 50 km radius 100 km off-center of the WD (cf. 
Fig. Ej left). 

While the spherical initial flame led to a maximum tempera- 
ture of 4.25 X 10^ K in the collision of the surface material and 
released 0.861 x 10^" erg prior to the collision, the irregular- 
shaped initial flame caused much more burning. It released 
1.14 X 10^*^ erg before the clash, and the maximum tempera- 
ture reached in the collision region was only 3.31 x 10^ K (cf. 
TableE). 

4.1.4. Irregular asymmetric ignition 



Due to the dipolar convection flow structure found in the 
pre-ignition phase by Kuhlen et al. (2006), a lop-sided flame 
ignition — ^possibly in many separate spots — is plausible. That 
ignition region may extend out to ~200 km and reach down to 
the center of the WD. Some overshooting to the opposite side 
may also be possible if the actual flow is multipolar The igni- 
tion kernels will quickly merge due to burning and the result- 
ing flame structure may look like a teardrop with an irregular 
surface. 

We examined two such configurations, one ignited on only 
one side of the WD (Model 2Tldl00, cf. Tabled, and another 
in which the ignition extended through the center (Model 
2T2dl00). These ignition configurations are shown in the 
right hand plot of Fig.0] 

While the one-sided ignition evolution proceeded in a 
way similar to the single-bubble ignition simulations (Model 
2B50d200c of Table [2 can serve as a reference simulation; 
its evolution is shown in Fig.[0, the extension of the ignition 
region to the opposite side of the WD's center had dramatic 
consequences. The burning did not ascend to the surface of 
the star on just one side, but evolved into two large irregular 
bubbles that moved in both directions (cf. Fig.|5}. Naturally, 
the collision of material sweeping around the surface occurred 
off-axis. In this set-up, two opposing effects altered the colli- 
sion strength. On the one hand, burning material on both sides 
of the star releases more energy and therefore the expansion 
of the star proceeds more rapidly. But, on the other hand, the 
collision occurs only slightly more than half way around the 
hemisphere, and therefore takes place earlier 

The reference Model 2B50d200c reached a collision tem- 
perature of 2.22 X 10^ K at a density of 1.41 x lO^gcm"^ 
and released 1.46 x 10^" erg in the burning. The one-sided 
teardrop ignition led to similar values — a maximum tempera- 
ture of 2.21 X 10** K at a density of 3.57 x 10^ g cm"^ and a 
nuclear energy release of 1 .25 x 10^° erg. 

For the two-sided teardrop ignition, considerably more en- 
ergy (2.30 X 10^*^ erg) was released in the burning prior to the 
collision. This decreased the collision strength dramatically. 
The temperature of the compressed fuel did not exceed 10^ K 
(cf. Table |3}, indicating that the effect of the earlier clash of 
the material coming from both poles cannot compensate for 
the enhanced expansion. However, it cannot be ruled out that 
multiple plumes breaking out at smaller angles with the cen- 
ter of the WD may collide more efficiently. But in this case 
the initial conditions have to be chosen carefully since bub- 
bles too close to each other will merge before reaching the 
surface. Looking at the evolution of our models ignited in a 
single bubble, this seems rather hard to achieve. The snapshot 
at f = 1 .0 s of Model 2B50d200c shown in fig.[T]indicates that 
by the time the burnt structure breaks our of the surface of the 
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TABLE 3 

Collision parameters for 2D simulations with dieferent ignition morphologies. Dashes mark cases where the density at which to measure the temperature is 

NOT REACHED IN THE COLLISION. 



Mode! 


resolution 


Tniax at coll. 


£nuc at coll. 


rmax(p>3xl0'5 


gcm-3) rn,„(p> lxlO^gcm-3) 


surface 






[10' K] 


[lO^^erg] 


at coll. [10' K] 


at coll. [10' K] 


detonation? 


2P50dl00 


384 X 768 


2.31 


1.14 


2.31 


0.507 


yes 


2B50dl00 


384 X 768 


4.25 


0.861 


4.25 


4.25 


yes 


2Tld200 


256x512 


2.21 


1.25 


0.991 




no 


2T2d200 


256x512 


0.958 


2.30 


0.0797 




no 




5.0-10" 1.0-10' 1.5-10= 2.0-10= 1-10= 2-10= 3-10' 4-10' 5-10= 2-10= 4-10= 0-10= 6-10= 1-10' 

X [cm] X [cm] x [cm] 

Fig. 5. — Evolution of an explosion simulation with the flame ignition extending to opposite sides of the WD (Model 2T2d200). 
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Fig. 6. — Maximum fuel temperature reached in the collision region as 
a function of the nuclear energy release prior to reaching that temperature 
according to Tables Q |2| (d ata points marked "2D"), Table |3] (data points 
marked "2Ds") and Tablel4l(data points marked "3D"). 

Star, it spans an opening angle of about 180° with the WD's 
center making the opposite-sided ignition the best choice. 



4.2. Lessons learned from 2D simulations 

Our exploration of setup parameters in 2D simulations re- 
veals that a key quantity determining the collision strength is 
the nuclear energy released on the flame bubble's way to the 
surface. Combining the data of Tablesmj] and|3] a clear cor- 
relation is visible between the maximum temperature reached 
in the collision of surface material and the amount of burning 
(cf. Fig.|5}. Whether accomplished through a change in res- 
olution, a displacement of ignition point, or a different mor- 
phology for the ignition region, less burning on the way out 
correlates with a stronger, hotter collision on the far side. This 
correlation arises naturally because the expansion of the star 
leads to the collision being spread out over a larger volume. 
The lower density also implies a greater heat capacity in the 
radiation field. On the other hand, more burning also implies 
more ash participating in the collision, which might make it 
stronger. But apparently the expansion effect dominates. 

A potential concern is that insufficient resolution of the 
compressed fuel region might make the temperature measure- 
ments unreliable. However, the clear correlation shown in 
Fig.|5]originates from simulations with different resolutions. 
The data points lining up well is an indication that the tem- 
perature measurement was credible even in the less resolved 
simulations 

The amount of fuel burned is a consequence of many uncer- 
tain aspects of the explosion physics, the specific algorithm 
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Fig. 7. — Evolution of a 3D explosion simulation with the flame ignition in a single bubble of radius 25 km displaced 100 km from the center of the WD (Model 
3B25dlOO). Top left panel: initial evolution of the flame front (blue to green isosurfaces coiTespond to ( = [0, 0.25, 0.35, 0.45] s). Other panels: later evolution 
with the logarithm of the density volume rendered and G = as blue isosurface indicating the flame front or, later, the approximate boundary between burned 
and unbumed material. 



used to implement the flame propagation, and the resolution 
of the simulation. For displacements that are not too extreme, 
stronger collisions are favored by increased distance of the ig- 
nition from the center. Alterations of the initial bubble shape 
diverging from the idealized spherical bubble model also have 
a substantial impact on the flame propagation. More complex 
initial flame shapes provide seeds for the developing nonlin- 
ear flame features. 

Two-sided ignition naturally burns more material, leading 
to greater expansion. This is only partially compensated by 
the earlier collision time. 



5. THE FULL STORY: 3D SIMULATIONS 

The entire WD was mapped onto a Cartesian computa- 
tional grid and, again, different ignition setups were tested. 
Compared with 2D, several general factors alter the collision 
strength in 3D. In the 2D-simulations, each flame feature cor- 
responds to a torus extending around the star Such com- 
plete burning does not occur in three dimensions, and so one 
might expect stronger collisions in 3D due to decreased ex- 
pansion. On the other hand, the additional degree of freedom 
also enhances the growth of the flame surface due to instabil- 
ities. This eff'ect has been seen in previous simulations where 
3D-, centrally-ignited setups released significantly more en- 
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Fig. 8.— Snapshot of Model 3B25dlOO at t = 3.3 s, as in Fig.Q but here 
with the temperature volume rendered instead of the density. 

ergy t han their 2D-counterparts jReinecke et alJ2002aHR6pk3 
|2005). This causes more expansion and weaker colHsions. 

Finally, when the collision geometry is no longer restricted 
by cylindrical geometry, one expects less focusing, again 
weakening the collision strength. 

5.1. Single-bubble ignition 

In two of the simulations, the flame was again ignited as 
a single spherical bubble. The subsequent evolution is given 
in Fig.0 Due to the interplay of burning and buoyancy, the 
burning bubble alters its shape from a sphere to a torus during 
the first few tenths of a second (cf., upper left pa nel of Fig.lTt. 
This evolution is very similar to the results of IZingale et alJ 
(|2005) who simulated a burning bubble on small scales fully 
resolving the flame structure. The developing torus is more 
regular than the respective flame structures found in the 2D 
simulations. Nonetheless, it is subject to instabilities, and 
irregular features grow, mostly on the outward side of the 
burning region. The seeds for these features are probably dis- 
cretization errors. 

Once the ashes reach the outer parts of the star (n.b., not 
necessarily the surface), they start to sweep around its core. 
Interestingly, the leading edge of this sweeping material is 
defined, even in 3D, by the former torus and therefore only 
slightly irregular (cf. upper right and lower left panels of 
Fig. 0. Consequently, the clash of the burned material still 
takes place in a well-defined spot on the opposite side of the 
WD (cf. lower right panel of Fig.0. In comparison with 2D 
simulations, this effect partially compensates for the lack of 
symmetry restrictions and makes the focus of the collision 
sharper than expected. In the collision, the temperature in- 
creases, as shown by the volume rendering of the temperature 
field in Fig.gl 

In the two simulations presented here, the initial flame bub- 
bles were displaced 100 km (Model 3B25dlOO) and 200 km 
(Model 3B25d200) from the center of the WD, respectively. 
Both simulations were carried out on a [512]^' cells compu- 
tational grid. To gain the maximum possible resolution, the 
fine-spaced inner uniform part of the computational grid was 



TABLE 4 

Collision parameters in the 3D simulations. 



Model Tniax at coll. at coll. p at coll. 
[10'' K] [10^° erg] [lO^igcm-^] 

3B25dlOO 0.892 2.79 < 1.3 x 10^ 

3P25dl00 1.34 1.01 < 1.5 x 10^ 

3P50dlO0 0.724 1.79 < 2.5 x 10^ 

3B25d200 no collision: WD unbound 

3Tld200 0.281 3.30 < 3.2 X lO'' 

3T2d200 no collision: WD unbound 



extended only slightly beyond the ignition radius, and there- 
fore the initial resolution of the flame was coarser in the model 
ignited further off-center Contrary to the 2D study where the 
initial resolution was kept constant while varying the ignition 
position, we had to sacrifice comparability between the sim- 
ulations to better resolution in the model ignited closer to the 
center. Because of the computational expense of 3D simula- 
tions, two parameters were changed at the same time, i.e., the 
distance of the ignition from the center and the perturbation 
imposed on the igniting bubble due to discretization errors. 

Unlike the 2D results, the flame ignited at 200 km off-center 
released more energy in 3D than the one ignited at 100 km 
off-center The reason, most likely, is the different discretiza- 
tion errors. The large features that develop at the end of 
the torus-dominated phase of the evolution were more pro- 
nounced in Model 3B25d200. Therefore more material was 
consumed and the energy release of 5.65 x 10^" erg was even 
sufficient to unbind the star (which has a binding energy of 
-5.20 X 10^*^ erg). In the simulation ignited 100 km off-center, 
the WD remained bound and the ashes breaking out of the sur- 
face swept around the core and collided on the opposite side 
as illustrated in Fig.0 

To gain more control over the initial perturbations that later 
affect the growth and buoyancy of the flame, two additional 
simulations were carried out where the ignition was, on the 
average, spherical, but actually composed of many smaller 
spheres of hot ash, so that the seeds for instabilities were 
present from the beginning. One, Model 3P25dlOO, had pa- 
rameters similar to Model 3B25dlOO, with the center of the 
ignition region positioned 100 km off-center, and the small 
flame kernels located inside a radius of 25 km. For the sec- 
ond. Model 3P50dlOO, the displacement of the center of the 
flame was unchanged, but the spherical aggregate of small 
bubbles filled a larger volume with radius of 50 km. In con- 
trast to the 2D simulation, the perturbations applied here de- 
creased the amount of burning (Table 1^. This can be under- 
stood from the the temporal evolution of the energy release 
and the bubble morphology. In the beginning, the energy re- 
lease in the highly perturbed case exceeds that for a smooth 
sphere, as expected from the faster development of its flame 
surface. However, this irregularity prevents the perturbed bub- 
ble from evolving a stable, toroidal structure. In terms of the 
overall energy release, the toroidal shape seems to be a favor- 
able configuration. For the spherical bubble models, as soon 
as irregular features form on top of the torus, the energy re- 
lease increases dramatically. This boost in the burning rate is 
much weaker for the perturbed case. 

In our 3D simulations, nuclear burning was suppressed 
three seconds after ignition. In all models but 3P50dlOO, it 
actually ceased earlier, since the fuel density in front of the 
flame fell below the threshold for burning, lO^gcm"^. For 
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Fig. 9. — Evolution of an explosion simulation with the flame ignition in a two-sided teardrop shape (Model 3T2d200). The first snapshot shows the initial 
flame (blue isosurface) with parts extending to both sides of the center of the WD. This leads to a two-sided evolution of the flame inside the WD, whose extend 
is indicated by the volume rendering of the logarithm of the density (center panel). The energy released in the thermonuclear burning is sufficient to overcome 
the gravitational binding of the star, giving rise to a weak explosion. In the right snapshot, the blue isosurface now indicates the approximate interface between 
burned and unburned material. This configuration is approaching homologous expansion and no collision between ash regions will take place at the surface of 
the star. 

Model 3P50dlOO, however, burning was still active at f = 3 s. 
Therefore, the nuclear energy release is only a lower limit for 
this calculation, and the maximum temperature in the colli- 
sion (found at f = 4.81 s) is an upper bound. 



TABLE 5 

Constraints on detonation initiations. 



5.2. Irregular asymmetric and dipolar ignitions 

As in the 2D simulations, irregular one- and two-sided 
teardrop-shaped initial flame setups were employed. The evo- 
lution of the two-sided flame ignition case (Model 3T2d200) 
is shown in Fig. O Once more, the flame propagates in two 
opposite directions. Unlike the 2D simulations, however, the 
burning releasing 8.16 x 10''"erg is suflicient and unbind the 
star The shape shown in the right panel of Fig. O marks the 
final stage of the evolution approaching homologous expan- 
sion. No strong collision of surface material occurred. 

Model 3Tld200, initiated with a one-sided teardrop-shaped 
flame, released about two-thirds of the WD's binding en- 
ergy. As with the other bubble ignitions, the burning material 
floated towards the surface, swept around the core of the WD, 
and clashed on the opposite side. The collision parameters are 
fisted in Tableg] 

5.3. Summary of the 3D simulations 

Tfie diversity of results found in 3D simulations is larger 
tfian tfiat in 2D simulations. Tfie measured quantities are sum- 
marized in Table 13 In two of tfie models tfie WD was even 
unbound. As plotted in Fig.|6l all peak collision temperatures 
in 3D simulations were lower tfian tfiose in 2D. 

A direct comparison between 2D and 3D models is diffi- 
cult since models witfi similar initial flames (in spfierical or 
teardrop-like sfiapes) release significantly more energy wfien 
performed in 3D. Tfie only way we found to lower tfie energy 
release fiere was by explicitly perturbing tfie initial bubble. 
Tfie two corresponding models, 3P25dlOO and 3P50dlOO, re- 
leasing 1.01 X 10^" erg and 1.79 x 10''°erg of energy in burn- 
ing, can be compared to 2D simulations releasing similar 
amounts of nuclear energy. Tfie closest examples would be 
2B50d200b (0.97 x lO^^erg) and 2B50dl00 (1.93 x lO^^erg). 
Tfiese acfiieved collision temperatures of 2.92 x 10^ K and 
2.00 X 10^ K, respectively, wfiile tfie collision temperatures 
for tfie 3D models were significantly lower (1 .34 x 10^ K and 
0.724 X 10''K). 



p[10*gcm"3] 


[10' K] 


M[g] 


R [km] 


detonation? 


10 


2.6 


2.5 X 10-3 


2 


no 


10 


2.7 


2.5 X 10-3 


2 


no 


10 


2.8 


2.5 X 10-3 


2 


yes 


10 


2.1 


2.0 X 10-5 


8 


no 


10 


2.2 


2.0 X 10^5 


8 


yes 


10 


1.8 


1.5 X 10^^ 


30 


no 


10 


1.9 


1.5 X 10" 


30 


yes 


3 


2.2 


2.0 X lO"** 


120 


no 


3 


2.3 


2.0 X lO-** 


120 


yes 


1 


2.4 


3.0 X lO" 


90 


no 


1 


3.0 


3.0 X 10" 


90 


no 


1 


3.0 


3.0 X 103" 


900 


no 



Two interpretations are possible fiere. Eitfier tfie collision 
temperatures are generally lower in 3D due to tfie additional 
degree of freedom decreasing tfie focusing, or tfiis may only 
be tfie case for models started witfi strong perturbations from 
a spfierical bubble since fiere tfie toroidal structure supporting 
focusing is suppressed. In tfie first case, tfie relation between 
collision temperature and energy release would be swallower 
for tfie 3D data points in Fig. |6l tfian tfie relation for tfie 2D 
data points. In tfie second case, tfie two 3D data points for 
tfie fiigfiest energy releases would fall onto tfie relation for tfie 
2D sample and tfie two low-energy 3D models would diverge 
from it. 

6. CONDITIONS FOR DETONATION 

6.1. Constraints on detonation ignitions in degenerate C+0 

matter 

In order to trigger a detonation a region must burn super- 
sonically and tfie size of tfiat region must b e larger tfian some 
critical mass ( Niemever & WooslevI (19971: Ipursi & Timmei 
I2OO61) . Because tfiat critical mass is very sensitive to tfie den- 
sity and composition, detonation becomes increasingly diffi- 
cult at low density and is sensitive to tfie carbon mass fraction 
in tfie unburned fuel. 

A series of calculations was carried out offline to study tfie 
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Fig. 10. — Maximum temperature reached in the collision region for fuel 
above the given density thresholds as a function of the nuclear energy release 
prior to reaching the maximum collision temperature (according to TableslTI 
121 and |3). The shaded and dashed regions con'espond to conditions where 
an initiation of a detonation is possible for temperatures reached in fuel of 
densities p > 1 X lO' gcm"^^ and p > 3 X 10^ gem"'', respectively. 

conditio ns for detonation usi ng the Kepler ID hydrodynam- 
ics code (Weaver et al. 1978). The procedure was identical to 
that described in Niemever & Wooslev ( 1997). 

A sphere composed of 50% by mass carbon and 50% oxy- 
gen of prescribed density, p, was given a temperature pro- 
file characterized by a central value, T^, and a linear decline 
over a specified range of mass, M, defining a radius, R, of the 
sphere. The sphere was then allowed to runaway inside of a 
much larger, cooler isothermal region to see if a successful 
detonation resulted. 

Table|5]gives the results. If the compression heats fuel with 
density above 10^ g cm"-' to a temperature over 1 .9 x 10^ K on 
a length scale of 10 km or more (the grid resolution) detona- 
tion will occur By 3 X 10^ g cm"-', the necessary temperature 
has risen to about 2.3 x 10 K on a scale of 100 km, and by 
p — 1 X 10^, it is impossible to detonate the star no mat- 
ter what temperature is achieved in the collision. The critical 
mass has become more than a substantial fraction of the en- 
tire star. Even though burning might occur with a supersonic 
phase velocity, detonation of carbon does not happen below 
10'' g cm"l 

6.2. Comparison with simulations 

Since it is the combination of temperature and density 
reached in a critical mass that decides whether detonation 
occurs, the maximum temperatures in the collision region 
were measured in our 2D simulations at a density exceed- 
ing 1 X 10' gem"-' and 3 x 10^ gem"-', respectively. Val- 
ues are given in Tables [2 El and |3] and the results are plot- 
ted in Fig. [TO] Densities in the compressed material above 
10 X lO'gcm"^ are only found in some rare cases in which 
the nuclear energy release prior to collision was lower than 
1.2 X 10^*' erg. For all simulations with an energy release 
lower than 2 x 10^" erg, the compression density of the un- 
burned material in the collision exceeded ~3 x 10^ g cm"^ and 
the temperature was determined there. 

From these measurements, we find that a detonation is ad- 
missible in some of the 2D models. In models with single- 
bubble ignitions displaced more than 200 km off-center, suf- 



ficient temperatures are reached at densities exceeding the 
threshold of 3 x lO^gcm"^. The spatial extent of the com- 
pressed region is not critical, since it is, at most, several grid 
cells in all cases (typical cell sizes in the collision region are 
several km). In models ignited at distances > 400 km off- 
center, temperatures above 2.5 x lO^K were reached even at 
densities above lO'gcm"^, rendering a detonation virtually 
certain. A bubble displacement around 200 km marks the bi- 
furcation value, where secondary parameters, such as resolu- 
tion, determine the feasibility of a detonation (cf. Table^. It 
seems unlikely that the flame ignition takes place at such large 
radii in SNe la, but similar results might come from a model 
with less efficient lateral burning ignited closer in. 

In all 3D simulations, the maximum compression tempera- 
tures in the collision (if one occurs at all) are too low to initiate 
a detonation. Moreover, these peak temperatures were found 
in material of densities falling short of the detonation thresh- 
old for at least one order of magnitude (cf. Tablel^. The com- 
pressed region did not reach densities above 10^ gem"-' near 
temperature maximum in any of the simulations, t he re asons 
thereof being the same as those discussed in Sect. 15.31 Even 
the 3D simulations releasing similar amounts of energy as 
some 2D simulations favoring a detonation did not reach the 
necessary densities. In the two-sided ignition Model 3T2d200 
the the effect of an increased energy release was even greater 
than its the 2D analog — it unbound the star. Thus, all 3D sim- 
ulations clearly fail to trigger a detonation. 

7. CONCLUSIONS 

The evolution of thermonuclear supernova models that ig- 
nite asymmetrically has been followed in two and three di- 
mensions. Parameters of the setup, such as the numerical res- 
olution, the displacement of the center of flame ignition from 
the center of the WD, and the ignition shape were explored in 
a systematic way in 2D, and the more interesting cases were 
explored in 3D. 

For all 2D simulations, the energy release in the nuclear 
burning falls short of unbinding the star. In the cases with one- 
sided ignition, the flame floats rapidly to the surface (see also, 
iLivne e t al.'2005). spreading laterally as it goes due to burning 
and instabilities. Since the star is still gravitationally bound, 
the emerging ashes sweep around the core of t he WD and 
collide on the opposite side, in agreement with iPlewa etall 
( 2004). The collision strength, and thus the maximum temper- 
ature reached in the compressed fuel correlate inversely with 
the nuclear energy released on the flame's way to the surface. 
The question of whether the compression of unburned mate- 
rial in the collision region is adequate to trigger a detonation 
was explored in detail, and the necessary criteria were set out. 
A detonation can only occur if the fuel temperature exceeds 
approximately 1.9 x 10^ K at a density above 10' gem"-'. At 
lower densities, detonation requires higher temperatures and 
eventually becomes impossible, for any temperature, for den- 
sities less than about 10^ g cm"-'. 

The conditions for initiating a detonation were met in sev- 
eral 2D calculations in which the flame ignited in a spheri- 
cal bubble more than 200 km off-center Less efficient pre- 
scriptions for the burning might have found similar conditions 
in simulations that ignited closer in. For the initial condi- 
tions and flame propagation model assumed, the results of 
Plewa et al. (2004), may be reasonable. However, since the 
flame model applied there is not based on a consistent treat- 
ment of the flame's interaction with the turbulent cascade, it 
is difficult to judge its validity. 
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In three dimensions, all simulations fell far short of initi- 
ating a detonation. Some even released sufficient energy to 
unbind the star There are several reasons for this difference. 
Lacking the artificial symmetry of 2D simulations, the focus- 
ing of the collision in 3D models can be weaker Indications 
for this were found in our simulations. At least two of the 3D 
simulations gave significantly lower collision temperatures 
than predicted by 2D models that burned similar amounts of 
fuel (see Fig.|5). A second, probably dominant effect is that 
3D models release more energy than their 2D analogs. This 
is mostly due to the additional degree of freedom in devel- 
oping flame surface area due to instabilities. Another contri- 
bution to the diff'erence may be the improved subgrid-scale 
turbulence model applied in the 3D simula tions. However , 
this eff'ect is expected to be minor ( Schm idt et a l] l2006bh . 
An asymmetrically ignited 3D model based on a different 
flame implementation (Calder et al. 2004) burned ~0.075 M© 
of material. This corresponds to a nuclear energy release of 
~1.2xl0^'^ erg — a value that falls in the range spanned by our 
parameter study. Therefore, although the Calder et al. model 
was not followed beyond the breakout of the ashes from the 
surface, one expects that that model would also have failed to 
trigger a detonation. 

While we found no example of a successful detonation 
in our 3D simulations, this possibility cannot be completely 
ruled out since the exploration of the parameter space was 
incomplete. An interesting possibility is that of double-sided 
ignitions, which may be possible in a teardrop-shaped ignition 
overshooting through the WD's center (or, in a simpler config- 
uration, as two opposed bubbles). Such a configuration short- 
ens the way the material has to travel towards the collision 
spot once ashes break out of the star's surface on both sides. 
Therefore the expansion of the star may not be as advanced 
as in the collision on the far side of a single-bubble breakout. 
In this case, higher temperatures and densities are expected 
in the compressed fuel. On the other hand, burning on both 
sides of the star releases more energy while the flames propa- 



gate towards the surface. This increases the expansion prior to 
collision. In our simulations, the latter effect dominated and 
the collision was weak in a 2D model. A similar 3D model be- 
came unbound and no detonation was found in either Thus, 
if a multiple surface breakout scenario is to work at all, no 
more than two widely separated ignition kernels are admissi- 
ble or there will be too much expansion. An open question is 
whether a special placement of two or three bubbles spanning 
a smaller angle than 180° with the WD's center might favor 
the first effect, increasing the collision strength by shorten- 
ing the path of the surface material. However, the bubbles 
cannot be too close or they would merge quickly due to burn- 
ing, without significantly compressing the unburned material 
between them, and there cannot be very many, or they will 
prematurely unbind the star 

Keeping in mind the uncertainties of the flame model and 
the incompleteness of the parameter space explored in 3D 
simulations, we conclude, that although a detonation due to 
the colliding surface material may, in principle, occur for 
certain — ^possibly artificial — ignition configurations, it cannot 
serve as a robust model for SNe la. The simulations presented 
here indicate that it may not be realized in nature at all. 

For models that remain gravitationally bound, failure 
to initiate a detonation will lead to pulsations of the 
WD star (fNomoto et al." '1976). This may be a second 
chance for triggering a detonation ( Arnett & Livne 19941; 
iBravo (feTTa rcfa-Senz 200^ and this occurrence will be ad- 
dressed in a follow-up study. 
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